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ABSTRACT 


Acoustic  wave  characteristics  in  the  upper  atmosphere  are  deduced 
from  temporal  and  spatial  electron  density  variations  obtained  with  the 
Arecibo  back-scatter  sounder.  A  theory  is  developed  which  relates  the 
observed  wave  characteristics  to  the  power  spectrum  of  the  wave 


sources. < 


ACOUSTIC  WAVES  IN  THE  IONOSPHERE 


I.  Introduction 

Observations^  of  the  lower  ionosphere  and  upper  D  region  (the 
region  from  80  km  to  150  km)  show  that  there  the  atmosphere  is  in 
a  state  of  random  motion.  Since  only  a  limited  amount  of  statistical 
information  can  be  obtained  frpm  rgg^gt  kerng  experiments,  it  is  natural 
to  attempt  to  use  ground  based  measurements  to  obtain  information  about 
the  random  state  of  thg  lc$y<3r  ionosphere.  This  VJQ?k  is  a  report  of  one 
attempt  to  obtain  a  reasonably  complete  physical  picture  of  one  class  of 
ionospheric  motions:  those  corresponding  to  random  internal  acoustic 
waves. 

The  radar  facility  at  Arecibo,  Puerto  Rigo,  offers  the  best  avail¬ 
able  spatial  and  temporal  resolution  of  the  backscattered  power,  and  the 
least  ambiguity  as  to  the  cross-section  of  the  scattering  process.  For 
altitudes  greater  than  110  km(  the  mean  free  path  of  the  gas  is  greater 
than  the  wave  length  of  the  radio  wave  (70  cm),  which  is  in  turn  greater 
than  the  Debye  length  of  the  gas  (10  cm).  Under  these  circumstances,  it 

can  by  shown  that  Thompson  scattering  occurs,  so  that  the  backscattered 

(2) 

power  is  proportional  to  the  electron  density  of  the  gas.  '  Furthermore, 

the  conditions  of  the  theory  are  approximately  met  down  to  about  80  km; 

electron  densities  determined  by  rocket  borne  measurements  are  in  fairly 

(3) 

good  agreement  with  the  radar  data.  '  Thus,  a  measure  of  electron 
density  is  available  in  the  80-150  km  region,  with  a  spatial  resolution  of 
about  6  km  in  the  horizontal  and  vertical,  a  time  resolution  of  10  *  sec. 
and  a  10  percent  accuracy  in  backscattered  power. 


We  shall  be  dealing  then  with  a  random  electron  density  which 
depends  on  spatial  position  and  time.  However,  a  further  simplification 
is  available  in  the  80-150  km  region,  when  data  are  taken  only  around 
noon.  It  is  thus:  if  the  phenomena  which  we  wish  to  observe  have  time 
scales  t  less  than  about  30  minutes,  we  may  suppose  that  the  electron 
density  variations  are  proportional  to  the  bulk  fluid  density.  For  example, 
to  have  diffusion  be  unimportant  in  a  cube  of  size  1  during  a  time  r,  we 

1 2  62 

must  have  But  a  typical  diffusion  coefficient  might  be  10  crrr/sec 

at  110  km,  if  1  is  6  km,  t  must  be  less  than  several  days.  Again,  if 

photo  chemical  reactions  are  present,  and  have  a  characteristic  reaction 

time  Tj^,  then  t.^>  t.  At  around  100  km,  might  be  one  second.  Thus, 

one  would  expect  the  electron  density  fluctuations  to  be  in  photochemical 

equilibrium.  When  the  sun  is  overhead,  this  equilibrium  implies  a 

constant  electron  concentration.  Finally,  it  has  been  shown  that  electron 

density  fluctuations  are  the  same  as  the  neutral  density  fluctuation  in  an 

(4) 

acoustic  wave  in  this  altitude  region'  ,  even  though  the  electron  velocity 
is  not  the  same  as  that  of  the  bulk  xluid,  due  to  magneto-ionic  effects. 

For  the  purpose  of  the  following  discussion,  we  assume  that  the 
backscattered  power  received  by  the  Arecibo  antenna  measures  the  bulk 
density  of  the  gas.  The  purpose  of  this  paper  is  to  describe  how  the 
statistics  of  the  bulk  density  fluctuations  can  be  used  to  determine  the 
motions  of  the  gas.  In  particular,  we  show  how  such  measurements  as 
described  above  can  be  used  to  deduce  the  source  strength  of  random 
acoustic  waves.  The  methodology  of  the  analysis  required  to  distinguish 
random  waves  from  other  kinds  of  random  motions,  such  as  turbulence, 
is  reviewed.  The  data  from  the  present  observations  is  presented.  A 
theory  is  presented  which  is  consistent  with  the  present  observations. 
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in  (3)  are  taken  along  T.  If  At  and  A3Tare  allowed  to  become  very  large, 

(51 

then  p  will  vanish  exponentially  rapidly'  ,  whereas  the  wave-like  com- 

*> 

ponent  of  R  varies  algebraicly 
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This  integral  may  be  approximated  by  the  method  of  stationary  phase, 
the  stationary  point  is 


grad  <0 
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Ale  =  Alt  , 
At  At 
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Provided  Ax]  At  are  sufficiently  large,  all  the  contributions  to  R  will 
come  from  wave  numbers  and  frequencies  which  satisfy  (4).  It  is  only 
when  Ax  and  At  become  large  that  the  wave-like  component  of  R  may  be 
separated  from  p  .  Hence,  we  obtain  the  following  rule:  If  the  antenna 
is  held  stationary,  the  waves  seen  will  be  those  vhose  group  velocity  is 
parallel  to  the  orientation  of  the  antenna,  s. 

Consider  now  a  still  atmosphere.  The  internal  waves  which 
may  propagate  in  such  an  atmosphere  have  the  following  general  proper¬ 
ties  (see  ref.  (7)  ):  For  s  oriented  at  an  angle  0  to  the  vertical,  there 
are  two  internal  gravity  waves  k^,  and  one  acoustic  wave  k^  which 

satisfy  equation  4.  Thus,  one  is  faced  with  determining  three  amplitudes 
from  R.  If  winds  of  unknown  magnitude  are  assumed  to  be  present,  then 
the  difficulties  are  magnified. 

However,  there  is  an  important  simplification  available  when 
the  antenna  points  directly  vertical.  In  this  case,  there  is  only  one 
solution  to  equation  (4),  which  is  an  acoustic  wave  propagating  vertically. 


This  is  the  case  studied  here. 

The  data  were  collected  in  December  1965  and  have  been  re¬ 
ported  in  ref.  (6).  The  two  main  features  of  the  data  which  we  wish  to 
describe  are  shown  in  Fig.  1,  which  shows  the  oscillatory  character  of 
the  correlation  function  for  large  Az,  At  and  Fig.  2.  Fig.  2  shows  the 
zeroes  of  the  correlation  function  for  large  lags.  The  error  bars  are 
estimated  as  in  ref.  (6). 

The  main  features  of  the  data  are:(a)  the  amplitude  of  the 
oscillations  does  not  increase  with  A  z,  and  (b)  the  zeroes  of  R  are  nearly 
straight  lines  parallel  to  the  Az  axis. 

III.  Theoretical  Model 

In  this  section  we  present  a  simple  model  of  vertically  propagating 
disturbance  which  has  the  following  properties.  1)  the  dispersion  re¬ 
lationship  yields  the  distribution  of  zeroes  shown  in  Fig.  2,  and  2)  the 
amplitude  is  nearly  constant  with  altitude  as  in  Fig.  1. 

Recall  that  in  a  constant  temperature  atmosphere,  an  internal 
wave  would  grow  exponentially  with  altitude,  provided  viscosity  were 
not  important.  It  is  known^  that  for  long  gravity  waves,  viscosity  is 
an  important  filter  in  the  100  km  range.  However,  such  waves  propagate 
nearly  horizontally,  and  hence,  at  a  given  altitude,  viscous  forces  have 
a  long  time  to  act  on  a  wave  packet.  Because  viscous  forces  at  a  given 
altitude  act  only  for  a  short  time  on  an  acoustic  wave  packet  propagating 
vertically,  it  seems  plausible  that  viscosity  would  become  important 
only  at  altitudes  higher  than  100  km. 

If  c  is  the  speed  of  sound,  p  is  fluid  density,  and  p  the  viscosity 
of  the  gas,  and  H  a  scale  height,  then  for  acoustic  waves,  an  appropriate 
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FIG.  2. 
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Reynolds  number  (the  ratio  of  viscous  to  inertial  forces)  is 

Re=  . 

V- 

Fig.  3  shows  that  this  parameter  is  quite  large  in  the  100-200  km  altitude 
range,  indicating  that  viscous  effects  are  unimportant  in  this  region. 

This  result  is  consistent  with  the  results  of  Maeda  (1964)  who  found  that 
acoustic  waves  were  damped  by  viscous  action  around  300  km  for  waves 
of  100  sec  period  (see  Fig.  11  of  Maeda  1964).  Since  the  period  of  the 
oscillations  in  R  is  nearer  300  sec,  Maeda’ s  calculation  suggests  that, 
these  disturbances  should  not  be  attenuated  by  viscosity  below  about 
400  km.  Furthermore,  the  detailed  calculations  of  Yanowitch  (1967) 
on  the  viscous  damping  of  the  acoustic  wave  may  be  used  to  compute  the 
expected  form  of  the  correlation  coefficient  R.  When  this  is  done,  one 
finds:  1)  the  damping  effect  switches  on  exponentially  fast  at  some  altitude, 
and  2)  above  that  altitude,  the  oscillations  in  R  increase  in  period. 

Neither  of  these  effects  is  observed. 

Fig.  5  shows  that  in  the  80-150  km  region,  there  is  a  strong 
temperature  gradient.  We  now  show  that  this  temperature  gradient  can 
be  invoked  to  explain  the  amplitude  dependence  of  the  data.  The  analysis 
below  is  analogous  to  that  given  by  Maeda  (1964),  except  that  low  frequency 
waves  are  considered.  If  D/Dt  is  the  substantial  derivative,  p  the 
pressure  and  u  the  fluid  velocity  taken  vertical,  then  the  equations  of 
motion  are 
Continuity 

1 

P 


R £ 

Dt 


9  u 
9  t 
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Momentum 


Du  , 

p  Dt  +  3  z 


+  g  =  0 


(2) 


Energy 


Dp  2  Dp 

Dt  ~  C  Dt 

By  linearizing  these  equations  and  cross  differentiating,  one  obtains 
(see  Maeda  1964) 

3  2u  _2  32u  8u 

- -  -  c  - -  +  Vg  -JT  =  0 


(3) 


(4) 


3t  3  z 

Here  c  ='\fyRTQ  ,  the  speed  of  sound,  g  is  the  acceleration  due  to 
gravity  and  Y  is  the  ratio  of  specific  heats, 

A  simple  approximate  solution,  valid  for  small  temperature 
gradients,  will  suffice.  Set 


u  =  A 


exp 


K 

2 


fZ  dz_ 

J  RT 
°  o; 


icot  A  /  v 

e  u  (z) 
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Then,  upon  substitution  into  equation  (4)  gives 
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i 
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o 
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We  assume  that  the  temperature  increases  exponentially  with  altitude 
(seeFig.4);  hence 


T  (z  =  80  km  )  exp  z/p 


-12- 


Using  this,  the  amplitude  factor  in  equation  (5), 


A  exp 


I  -gP 

I  2RT(0) 


\ 

exp  -  z  /  p  j  , 
/ 


(7) 


varies  by  at  most  50%,  as  1<  z/P  <  2  in  the  120-150  km  altitude  range. 

A 

Here  T(0)  =  Tq(z  =  80  km).  In  this  range,  u  has  the  approximate  form 

u  ~  e^Z ,  k  =  ~\l  J  1  1  (8) 

VvRTc  -  4H2  -2HF 


Here,  H,  the  scale  height,  is  RTQ/g.  As  P  ~  40  km,  and  H  ~  8  km,  it 
is  seen  that  the  main  effect  of  the  temperature  gradient  is  to  reduce  the 
exponential  increase  in  amplitude  of  the  wave  while  leaving  the  dispersion 
relationship  nearly  unchanged. 


IV.  .The  Power  Spectrum 

In  this  section  we  use  the  approximate  solution  of  the  preceding 
section  to  estimate  the  form  of  the  correlation  functions  of  the  waves. 
According  to  equation  (7),  (8),  the  disturbance  in  density  (which  is 
proportional  to  u'due  to  linearity)  has  the  form 
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^  /  ft 
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-CO 


exp 


z  * 

Pi 


icot  -ikz  . 
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where  k  is  given  by  equation  (8).  Forming  the  correlation  coefficient 
R,  we  have 

„  0  °°  ^  ,  \ja*/  -  *  .  iwAt  -ikAz 

R  =  \  <  dA  (o,  z)  dA  (w,  z  +  Az)  >  e  e 


exp 
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2RT 
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1  z  \  , 
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In  the  limit  of  large  At,  Az  this  integral  has  a  stationary  point  at 


d  k  _  A  t 
d  co  **  A  z 

which  can  be  shown  to  be 

A  t  _  to 

A  z-  v  RT  k 

T  o 

In  terms  of  frequencies  this  result  may  be  written  as 


2  2 

Referring  now  to  figure  3,  it  is  clear  that  (yRTy)  (At)  >>  (Az)  , 
so  that  we  may  expand  the  formula  for  co  as 


According  to  this  model,  the  acoustic  cutoff  frequency  in  a  temperature 
gradient  is 


which  is  lower  by  a  factory(l  p-  )  from  the  isothermal  case.  Using 
formula  (9)  one  can  associate  oscillations  occurring  in  R  at  given  values 


The  amplitude  of  the  oscillations  is  easily  worked  out  to  be 


R  =  P(o),  z,  Az)  A  (z)  cos 


f 

(At)2  yRTo  -  (Az)2 

n  2H  i 
(!  — ) 

1/2  i 

2 

4ET 

+  4  fy 

/ 

where  A  (z)  has  the  form 


6Xp  2RT  exp  "  f  ] 
o 


For  (Atr  yRT q> > Az , the  amplitude  of  the  oscillations  may  be  approxi¬ 
mated  by 


R  =  A(z)  P(co,  z,  Az)  cos/^Y*,RTo  (1  ~'T~)  +  ?  V  {11) 


It  will  be  seen  shortly  that  the  observed  oscillations  are  very 
close  to  the  acoustic  cutoff  frequency.  In  this  case,  the  lines  of  constant 
phase  are  straight  lines,  parallel  to  the  Az  axis  .  The  period  of  these 
oscillations,  taking  account  of  the  finite  temperature  gradient,  is  about 
310  sec  at  100  km.  The  observed  zeroes  of  R  (lines  of  constant  phase) 
(Fig.  2)  form  straight  lines  parallel  to  the  Az  axis,  and  their  mean 
spacing  yields  an  observed  period  of  320  sec. 

To  construct  a  power  spectrum  of  the  source  strength  of  the 
waves,  a  value  of  Az,  At  is  selected,  and  the  amplitude  of  this  is 
estimated  from  the  data.  This  amplitude  is  a  measure  of  the  power 
spectrum  P(co,  z,  Az).  The  frequency  is  found  by  equation  (10), 


rwraias 


Using  this  method,  the  results  shown  in  Fig.  6  have  been 
deduced.  It  is  reassuring  to  note  that  the  amplitudes  for  various  A  z, 

A  t,  reduce  to  a  single  curve  when  plotted  this  way. 

V.  Discussion 

One  might  compare  these  results  with  a  simple  nonpropagating 
disturbance,  of  the  form  u~e1C0-^^.  Experimentally,  the  correlation 
function,  R,  should  behave  as  independent  of  z.  Actually,  in  the 

region  of  A  z  ,  At  space  where  the  oscillations  are  observed,  this  is 
nearly  the  case.  Clearly,  the  equations  which  would  describe  small 
oscillations  and  the  acoustic  cutoff  frequency  are  the  same  as  those 
presented  in  Section  III.  For  all  practical  purposes,  then,  these  waves 
are  very  nearly  standing  waves. 

Because  of  this  fact,  a  calculation  of  the  energy  flux  carried  by 
the  wave  is  uncertain,  involving  the  difference  of  two  small  numbers. 

The  present  estimate  is  1.4x  10  watts/rm.  This  value  is  comparable 
with  that  predicted  by  Hines  (1965)  for  the  energy  flux  due  to  internal 
gravity  waves.  Although  it  can  be  shown  that  the  disturbance  from  a  point 
source  of  energy  in  the  lower  atmosphere,  ignoring  the  ground  effect, 
is  equally  divided  between  acoustic  and  gravity  waves,  it  does  not 
necessarily  follow  that  the  energy  flux  at  100  km  from  these  two  wave 
modes  is  equal.  The  present  data  do  suggest  this. 

Another  explanation  of  the  data  supposes  that  the  disturbances 
observed  have  their  origin  in  internal  gravity  waves  at  a  ground  level. 
Consider  an  internal  wave  very  near  the  gravity  wave  cutoff  frequency 
at  ground  level,  propagating  upward.  As  such  a  wave  encounters  the 
temperature  variations  in  the  atmosphere,  it  will,  in  some  regions, 


become  damped.  We  suppose  that  as  the  wave  propagates  through  the 
temperature  field  of  the  atmosphere,  a  mode  conversion  process  occurs, 
in  which  the  damped  gravity  wave  becomes  an  acoustic  wave  when  the 
temperature  drops  so  far  as  to  allow  the  disturbance  to  propagate  as  an 
acoustic  wave. 

According  to  this  model,  the  power  spectrum  of  the  waves 
should  be  bounded  above  by  w  q  «  the  gravity  wave  cutoff  frequency  at 
ground  level,  and  below  by  o  .  the  acoustic  cutoff  frequency  at  100  km. 
A  calculation  on  these  lines  shows  that  one  can  adequately  explain  the 
width  of  the  spectrum  observed  with  these  ideas.  However,  one  would 
expect  that  the  energy  flux  would  be  an  order  of  magnitude  down  from  that 
of  internal  gravity  waves  with  this  model.  It  is  not  clear  if  this  is  the 
case. 

Finally,  it  should  be  emphasized  that  the  methods  presented 
here  may  be  used  to  analyze  internal  gravity  waves  also;  provided  the 
antenna  is  directed  at  an  angle  to  the  vertical. 
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